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FOREWORD 


The  state  of  the  earth's  upper  atmosphere  under  normal  and 
disturbed  conditions  is  of  considerable  importance  for  a  variety 
of  military  applications,  ranging  from  ionospheric  radio  wave  pro¬ 
pagation  to  high-altitude  nuclear  explosions.  The  biggest  reservoir 
of  energy  storage  in  the  atmospheric  E-region  is  provided  by  the 
Schumann- Runge  photodissociation  of  molecular  oxygen,  Which  yields 
one  oxygen  atom  in  the  excited  ^D  state.  This  excited  atom  is 
rapidly  de-excited  ("quenched")  by  collision  with  atmospheric  molec¬ 
ular  nitrogen.  The  question  of  the  possible  vibrational  excitation 
of  nitrogen  by  this  process  and  of  its  dt -excitation  by  tne  vibration- 
translation  energy  transfer  between  N2  and  0  (vrtiich  is  also  studied 
here)  is  a  very  important  one,  since  the  atmospheric  E- region  is  far 
from  thermodynamic  equilibrium,  and  the  effective  vibrational  popula¬ 
tion,  vhich  is  partly  determined  by  these  two  processes,  may  affect 
both  the  neutral  and  ionic  atmospheric  chemistry  (in  particular  the 
de-ionization)  as  well  as  the  airglow. 

lhe  present  calculation  of  these  reactions  provides  the  basis 
for  an  experimental  investigation  of  these  same  processes  vhich  is 
now  being  conducted  at  Stanford  Research  Institute  under  ARPA  contract; 
the  present  calculation  guides  and  complements  those  measurements  as 
well  as  providing  a  basis  for  further  experimental  and  theoretical 
work  in  this  area. 

This  work  was  performed  under  sponsorship  of  the  Strategic 
Technology  Office  of  the  Advanced  Research  Projects  Agency. 
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ABSTRACT 


As  an  extension  of  earlier  work  on  the  quenching  of  electroni¬ 
cally  excited  alkali  metal  atoms  by  N2  and  CO,  which  uses  a  model  of 
curve  crossing  involving  an  ionic  intermediate  state,  here  we  treat 
the  quenching  of  0(^D)  by  N2,  a  process  which  proceeds  by  curve 
crossing  of  covalent  states.  The  model  is  also  applied  to  the  vibra¬ 
tional  relaxation  of  N2  in  collision  with  0(^P)  atoms  and  to  the  uni- 
molecular  decomposition  of  N20.  The  order  of  magnitude  of  the  rate 
constants  for  all  three  processes  can  be  explained  by  the  present 
model,  using  the  same  value  of  the  electronic  coupling  potential  for 
the  curve-crossing  process;  this  value  is  significantly  larger  than 
the  corresponding  value  for  isolated  0(3P)  and  O^D)  atoms.  In  the 
quenching  of  OC^D)  by  N2,  the  electronic-vibrational  coupling  is  weak, 
channeling  less  than  5  percent  of  the  electronic  quantum  into  vibra¬ 
tional  levels  of  N2,  but  the  quenchant  N2  is  likely  to  be  rotation- 
ally  excited.  The  vibrational  relaxation  of  N2  by  oxygen  atoms 
arising  due  to  curve  crossing  rather  than  to  adiabatic  ("Landau- 
Teller")  collisions  has  a  significant  activation  energy,  ~  0.8  ev,  so 
that  the  rate  coefficient  for  vibrational  relaxation  due  to  curve 
crossing  is  expected  to  show  a  strong  temperature  dependence  and  to 
dominate  at  temperatures  above  600°K. 
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I.  INTRODUCTION 

In  collisions  of  atoms  and  molecules  at  elevated  temperatures, 
the  flow  of  energy  between  electronic  and  vibrational,  rotational, 
and  translational  states,  i.e.,  quenching  processes,  can  have  signifi¬ 
cant  effects  on  both  transient  and  steady- state  molecular  properties. 
Recently,  calculations  were  presented  on  the  quenching  of  electronically 
excited  alkali  atoms  by  simple  diatomics  (Refs.  1  and  2),  e.g., 

Na(32P)  +  N2  -♦  [Na+  -  N“)  -♦  Na(32S)  +  N2(v)  +  K.E.  (1) 

vdnere  it  was  shown  that  the  vibration-electronic  coupling  was  strong 
due  to  the  unstable  intermediate  state,  vdiich  is  largely  ionic  in 
character.  As  an  outgrowth  of  this  work,  it  was  recognized  that  a 
similar  curve-crossing  mechanism  would  be  operative  in  the  relaxa¬ 
tion  of  vibrational  energy  in  N2  during  collisions  with  alkali 
atoms; 

N2(v)  +  Na(32S)  -♦  {N2  -  Na+)  -♦  N2(v '  <  v)  +  Na(32S)  +  K.E.  (2) 

These  calculations  (Ref.  3)  demonstrated  that  this  vibrational 
relaxation  process  may  have  significantly  larger  cross  sections  than 
other  vibration-translation  energy  transfer  processes  vdiich  do  not 
involve  a  curve  crossing. 

Reactions  (1)  and  (2)  proceed  by  the  crossing  of  potential 
energy  curves  (a  Landau-Zener  process)  through  an  ionic  intermediate 
state  (Ref.  4),  although  the  quenching  model  developed  in  the  earlier 
studies  is  not  restricted  to  the  crossing  of  ionic  and  covalent  states. 
In  this  paper,  a  quenching  reaction  is  investigated  vMch  proceeds 
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entirely  through  the  crossing  of  covalent  potential  energy  curves* 

By  analogy  to  the  ionic  quenching  studies,  the  following  reactions 
all  proceed  through  the  same  or  similar  potential  energy  curve 
crossings  (Ref.  5); 

N20(X1£'f )  -*  NjCX^E*)  ♦  0(3P)  (unimolecular  decomposition)  (3) 
0(-4»  +  N2  -*  {N20}  -*  0(3P)  ♦  H2(v,  J)  +  X.E.  (quenching  of 

0(1D>  by  H2)  (4) 

0(3P)  +  N2(v)  -*  (N20)  -♦  0(3P)  +  N2(v'  <  v)  K.E. 

(N2  vibrational  relaxation  by  0  atoms).  (5) 

Thus,  in  the  present  work  a  simple  but  unified  analysis  of 
reactions  (3)  -  (5)  is  presented, based  on  the  same  electronic  states 
of  N20.  In  Section  II  an  analysis  is  given  of  the  potential  energy 
curves  of  low-lying  electronic  states  of  the  NjO  system,  and  then  in 
Section  HI  the  nature  of  the  quantum  mechanical  collision  problem, 
in  particular  the  curve  crossings  and  the  magnitude  of  the  transi¬ 
tion  matrix  element,  is  discussed.  The  specific  processes  (3),  (4), 
and  (5)  are  discussed  in  Sections  IV,  V,  and  VI,  respectively,  and  a 
concluding  discussion  is  given  in  Section  VII. 

The  present  covalent  problem  has  considerable  differences  from 
the  previous  ionic  problem,  arising  from  the  differences  between 
short-range  (valence)  versus  long-range  (Coulomb)  interactions.  For 
an  ionic  intermediate,  a  long-range  crossing  (Ref.  6)  occurs  due  to 
the  Coulomb  interaction,  e  /R;  and  the  transition  matrix  element  at 
the  crossing  is  a  function  of  the  crossing  distance  (Refs.  7  and  8). 
Further,  because  of  the  long-range  interaction,  the  distance 
R(Na-N2)  or  R(Na+  -  Nj)  is  significantly  greater  than  the  distance 
R(N  -  N),  and  thus  the  geometrical  orientation  of  the  Na  relative  to 
N2,  or  the  angle  N-ft-Na,  is  relatively  unimportant.  However,  for 
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the  covalent  problem  R(N„  -  0)  =  R(N-N),  and  thus  the  orientation, 

*  A 

as  represented  by  the  angle  N-N-O,  is  of  major  importance  for  the 
potential  energy  surfaces.  Further,  the  curve  splitting  or  transition 
matrix  element  will  presumably  depend  on  effective  quantum  numbers  and 
molecular  symmetries  rather  than  on  the  crossing  distance.  Also, 
because  of  these  orientational  factors,  rotational  excitation  may  be 
particularly  important.  These  general  characteristics  of  ionic 
versus  covalent  crossings  are  of  importance  if  one  wishes  to  extend 
the  analysis  to  systems  other  than  NgO. 

It  is  also  appropriate  to  note  the  philosophy  of  the  present 
work,  vrtiich  is  to  identify  problem  areas  in  a  given  class  of  important 
reactions  by  a  semi- empirical  approach,  rather  than  to  give  a  quanti¬ 
tative  treatment.  Such  quantitative  treatments  will  become  possible 
within  the  next  several  years,  with  the  development  of  powerful 
numerical  techniques  (Ref.  9),  but  at  present  there  is  still  a  signifi¬ 
cant  gap  between  experimental  data  and  understanding,  vhich  can  be 
usefully  filled  by  qualitative  models  such  as  that  developed  in  this 
work. 

In  addition,  it  is  important  to  note  that  reactions  (4)  and  (5) 
are  of  considerable  aeronomic  importance.  The  quenching  of  0(^D)  by 
atmospheric  N2  is  a  potentially  important  source  of  vibrational  energy 
in  the  upper  atmosphere  (Refs.  10  and  11),  vhile  the  relaxation  of 
N2  vibrational  energy  by  0  atoms  may  be  an  important  loss  mechanism. 
Measurements  are  available  on  both  these  reactions  vhich  show  that 
the  quenching  rate  coefficients  are  large  (Refs.  12  and  13),  of  order 
1/4  gas  kinetic,  while  0  atoms  are  known  to  be  effective  in  relaxing 
N2  vibrational  energy  (Ref.  14)  at  least  in  high- temperature  systems. 

The  following  conclusions  can  be  drawn  from  the  model  developed 
in  this  paper:  the  vibration-electronic  coupling  in  the  quenching 
of  (X^D)  by  N2  is  found  to  be  weak,  channeling  only  about  5  percent 
of  the  electronic  quantum  into  vibrational  levels  of  the  quenchant 
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N2,  the  quench  ant  N2  is  likely  to  be  rotationally  excited  due  to 
the  importance  of  noncollinear  configurations,  and  the  relaxation 
of  N2  vibrational  energy  by  0(3P)  atoms  is  accompanied  by  a  signifi¬ 
cant  activation  energy,  ~  0.8  ev,  which  indicates  a  strong  tempera¬ 
ture  dependence  in  extrapolating  measured  high-temperature  rate 
coefficients  (Ref.  14)  to  lower  temperatures. 


II.  POTENTIAL  ENERGY  CURVES 

The  ground  state  of  the  triatomic  molecule  N20  is  known  to  be 

a  linear  structure  (N-N-O)  of  symmetry  ^T+,  so  that  it  clearly  dis- 
.  1  +  1 

sociates  into  N2(X  +  0(  D)  in  order  to  conserve  spin.  On  the 

other  hand,  O(T))  lies  1.98  ev  above  0(3P),  so  that  at  infinite 

separation  the  system  N0(X^E*)  +  0(3P)  lies  1.98  ev  below 
1  +  1  *  y 

N2(X  Eg)  +  0(  D).  There  is  thus  no  question  that  a  curve  crossing 
does  occur.  In  fact  the  situation  is  rather  complex:  Peyerimhoff  & 
Buenker  (Ref.  16)  have  made  quantum  mechanical  configuration  inter¬ 
action  calculations  for  the  system  N-N-0  as  a  function  of  the  angle 
a  =  N-N-O,  at  the  equilibrium  internuclear  separations  rg(N-N)  = 

1.128  A  and  rg(N-0)  =  1.184  A.  It  should  be  noted  that  the  state 
assignments  for  electronically  excited  states  of  N20  as  listed  in 
Herzberg  (Ref.  15)  do  not  agree  with  the  calculations  of  Peyerimhoff 
and  Buenker  (Ref.  16);  since  the  assignments  in  Herzberg  are  based  on 
very  limited  data,  it  appears  preferable  to  use  the  analysis  of 
Peyerimhoff  and  Buenker. 

For  a  collinear  configuration,  a  =  N-ft-0  =  180°,  as  the  N2(X^Eg) 

and  0(3P)  systems  are  brought  together,  there  are  two  electronic  states 

3E~  and  3n  of  successively  higher  energy;  when  N„(X1E+)  and  O('h))  are 

1  +  l  ^  l  " 

brought  together,  there  are  three  states  X  £  ,  a»  n*  All  these 
states  are  shown  in  Fig.  1  for  R(N-N)  =  rg(N-N)  =  1.128  A.  It  should 
be  noted  that  except  for  the  X*£+  state,  the  curves  for  the  other 
states  have  been  faired  in  using  the  energy  value  from  Ref.  16  for 
R(N-O)  =  re(N-0)  =  1.184  A. 

At  small  internuclear  distances,  the  and  3n  states  become 
Rydberg,  i.e.,  highly  repulsive,  while  the  and  3£"  states  are  of 
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POTENTIAL  ENERGY  («v) 


FIGURE  I .  Potential  Energy  Cun*  for  Um- Lying  $tat«  <rf  H»  N,-0 
System.  Co  Mi  near  Configuration  2 
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valence  character  and  thus  somewhat  less  repulsive.  The  X^E+  ground 
state  of’NjO  has  been  characterized  by  a  Morse  potential  using  numeri¬ 
cal  values  from  Herzberg  (Ref.  15).  The  activation  energy  associated 
with  the  crossing  of  the  1£+  -  ^E curves  has  been  taken  from  the  high- 
pressure  limit  of  the  unimolecular  decomposition  of  N20(X^E+)  into 
0(3P)  +  N2(X1Eg) — see  Ref.  18.  The  details  of  this  are  presented  in 
Section  IV,  together  with  a  brief  review  of  pertinent  parts  of  the 
decomposition  reaction.  Although  it  is  not  possible  to  give  a  relia¬ 
ble  estimate  of  the  accuracy  of  most  of  these  potential  energy  curves, 
it  is  unlikely  that  errors  in  relevant,  portions  of  the  crossing  region 
exceed  1  ev,  at  least  for  the  collinear  configuration,  a  =  180°;  the 
decrease  in  the  well  depth  of  the  ground  state  of  N20  with  a  decreas¬ 
ing  from  180°  is  taken  from  Peyerimhoff  and  Buenker  (Ref.  16). 

For  a  noncollinear  configuration,  a  twofold  degeneracy  existing 

for  states  other  than  E  states  is  removed,  so  that  each  n>  A . 

state  splits  into  two  states  of  lower  symmetry:  for  the  appropriate 
point  group  C  these  are  described  as  A'  and  A',  respectively--see 
Herzberg  (Ref.  15),  pp.  10,  576.  In  general,  the  energies  of  these 
different  states  will  vary  differently  with  a,  with  R(N-N),  and  with 
R(N-O);  a  possible  schematic  diagram  is  shown  in  Fig.  2  (Ref.  19). 

It  must  be  stressed  that  the  details  of  Fig.  2  are  certainly  not  cor¬ 
rect,  but  it  does  show  the  large  number  of  possible  curve  crossings 
whose  number  and  position  can  vary  quite  significantly  with  cr,  R(N-N), 
and  R(N-O) . 
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FIGURE  2.  Potential  Energy  Curves  for  Low-Lying  States  of  the  N^O  System: 

Bent  Configuration,  Corresponding  toflr  =  N-ft-0  =  120°.  State 
Assignments  are  in  Parentheses;  P2),  etc..  Correspond  to  the 
Linear  Configuration,  a  =  180°,  of  Fig.  I . 
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III.  THE  CURVE-CROSSING  MODEL 


A.  THE  COORDINATE  SYSTEM  AND  THE  HAMILTONIAN 

Three  different  modes  of  motion  are  considered  here:  namely, 
electronic,  vibrational,  and  translational/rotational.  Iftese  will  be 
described  in  terms  of  the  coordinates  r,  £,  and  R,  respectively, 
defined  as  shown  in  Fig.  3  (for  clarity  only  one  electron  is  shown). 
This  definition  has  the  advantage  that  the  motion  of  the  center  of 
mass  is  eliminated,  and  the  kinetic  energy  terms  in  the  Hamiltonian 
are  separable  and  are  given  by  the  sums  of  squares  of  the  appropriate 
momenta.  Thus, 


H(o>  -  T(r)  +  T<£)  ♦  T(R)  +  V(r,  £,  R)  (6) 

2  2 

vftere  T(<j)  =  -(ft  /2M)v^  are  the  appropriate  kinetic  energy  operators 
in  terms  of  the  proper  reduced  mass  M  and  the  coordinate  £  (=  r,  £,  R), 
and  the  potential  energy  function  V(r,  £,  R)  involves  coupling  of  all 
three  variables.  If  the  following  ansatz  for  the  wave  function  is 
made, 

^  <£!  £>  Vs)  \  <W  5.)  fx<£)  •  <7> 

vftere  cp.  (r;  R)  is  the  spatial  part  of  the  electronic  wave  function, 
vftich  depends  parametrically  on  the  coordinate  R  (dependence  on  the 
vibrational  coordinate  £  is  ignored);  vrtiile  (s)  is  the  appropriate 
spin  wave  function;  8,(p;  R)  is  the  vibrational  wave  function  which 
again  depends  parametrically  on  R;  and  f  (R)  is  the  translational- 

— '  A  ~ 

rotational  wave  function,  describing  the  motion  of  0  relative  to  Ng. 
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FIGURE  3  .  The  Coordinate  System 


The  curve  crossings  all  involve  a  change  in  spin  quantum  number, 
so  that  they  arise  from  the  spin-orbit  coupling  term  in  the  Hamiltonian, 

H50  *  5(r)  (L.s)  (8a) 

S(r)  =  (l/2m2c2)(l/r)(dV/dr)  -  (Ze2/2m2c2)(l/r3)  ,  (8b) 

vhere  L  is  the  angular  momentum  operator  of  the  electron  and  s  is 
its  spin  operator,  and  V(r)  Ze2/r  is  the  central  field  potential 
energy. 

B.  LANDAU- ZENER  FORMULA  AND  TRANSITION  MATRIX  ELEMENT 

Following  the  characterization  of  the  potential  energy  curve  in 
Section  II,  the  probability  of  a  diubatic  curve  crossing  must  be  cal¬ 
culated.  For  two  curves  corresponding  to  states  \>  \i  vhich  intersect 
at  R  =  5c>  this  probability  is  assumed  to  be  given  by  the  Landau-Zener 
expression  (Ref.  20), 
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(9) 


expl'VRc)}’ 


vrfiere 


W 


iy<R>'2 

lw«MV  -  VI 


(10) 


1  /2 

with  u  =  {2K(R  )/M]  '  -  relative  velocity  at  the  crossing  point 

c  c 

Rc  at  vdiich  the  relative  kinetic  energy  is  K(Rc>. 

To  determine  the  transition  matrix  element  at  the  crossing 
point,  V^(RC)»  note  that  the  transition  is  essentially  an  electronic 
one;  however,  there  is  the  possibility  of  a  change  in  intemuclear 
separation  p  of  the  N2  molecule  or  more  correctly  of  a  change  in 
effective  .vibrational  quantum  number  of  the  N-N-0  system  vdiich  corre¬ 
sponds  asymptotically  to  a  change  in  the  vibrational  quantum  number 
of  N2. 


Thus,  with  the  ansatz  (8)  for  the  spin-orbit  Hamiltonian  and 

the  ansatz  (7)  for  the  total  wave  function  of  the  system  (from 

vhich  the  translational  wave  function  f,(R)  has  to  be  taken  out  in 

X  ~ 

order  to  arrive  at  the  transition  matrix  element  V^(RC)  t0  be 
inserted  in  the  Landau-Zener  formula),  the  following  equation  is 
found 


vv = 


,so 


cp,(r;  R)X  (s)9.  ( p;  R),  How  cp(r;  R)Xm  (s)0  ,(p;  R) 
~  ~  X  X  *  ~  (i  ~  ~  (J  pi  ~ 

1/2 


=  yeVv 


WRc>  ■  • 


R  =  R, 


(11) 
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where 


Vel,  (R  )  = 
c' 


SO 


R)Xx(s)  ,  H  tpu(r;  R)Xw(s) 


«  Xu<R=>  ' 


9X<£!  £>>  6L!<£5E> 


R  =  R. 


R  =  R. 


(12a) 


(12b) 


In  other  words,  the  transition  matrix  element  is  the  product  of  an 

’el 

electronic  transition  matrix  element  V  ^  multiplied  by  a  vibra¬ 
tional  overlap  integral  (written  as  the  square  root  of  a  "collision- 
modified  Franck-Condon  factor"),  both  evaluated  at  the  crossing 
point. 

To  lowest  order, 

yel^(Rc)  -  Vel^(R  *♦  •)  =  Vel  {OC^  -  3P2)}  (13a) 

for  the  electronic  transition  matrix  element,  and  for  the  isolated 
oxygen  alone,  Yamanouchi  and  Horie  (Ref.  21)  give  the  result, 

V®1  0(^2  -  2P2)  -  80  cm-1  .  (13b) 

Accordingly,  one  would  expect  the  transition  matrix  element  between 

the  ^£+  and  triplet  curves  to  be  of  this  order  (Ref.  19).  To  obtain 

detailed  numerical  values  of  the  transition  matrix  elements  would  be 

a  major  computational  problem,  and  thus  our  approach  has  been  to  use 

various  experimental  data  to  obtain  plausible  empirical  values  for 
el 

V  ,  which  is  treated  as  a  parameter.  This  point  is  discussed  further 
in  Sections  IV  and  VII. 

C.  THE  COLLISION-MODIFIED  FRANCK-CONDON  FACTOR  qC 

To  the  approximation  of  Eq.  (11),  the  effect  of  the  vibration  on 
the  transition  matrix  element  enters  solely  through  the  vibrational 
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overlap  integral  and  is  customarily  expressed  in  terms  of  its  mean 

r 

square,  the  collision-modified  Franck-Condon  factor  q  (R)  of  Eq.  (12b). 


In  the  limit  of  infinite  intemuclear  separation  R, 
just  the  ordinary  Franck-Condon  factor  q^ 


Xu 


q~  (R)  becomes 
Xu 


q  .  ( R  -*  ®)  :■  q, 

H  Xu  MXU 


eu(£> 


(14) 


Here  0^(£)  is  the  vibrational  wave  function  for  an  isolated  N2  mole¬ 
cule  and  thus  the  Franck-Condon  factor  q  is  just  a  Kronecker  fl; 

Af* 


1  if  v  =  v 
X  U 

0  if  VX  *  VM 


(15) 


where  v^,  v^  are  the  vibrational  quantum  of  N2  in  states  \, 

To  obtain  an  estimate  for  the  possible  vibrational  excitation  or 
de-excitation  in  curve-crossing  collisions  using  the  present  model, 
we  have  to  go  beyond  the  approximation  of  q  ,  and  in  fact  the  pre¬ 
vious  work  of  Nikitin  (Ref.  22)  and  Bauer  (Ref.  23)  is  extended  here. 


Because  of  the  considerable  difference  in  the  interaction  of 
N2  with  0(1D)  and  0(3P),  as  shown  in  Figs.  1  and  2,  one  would  expect 
that  the  collisional  distortion  of  the  vibrational  motion  will  be 
different  for  the  two  states,  (which  we  will  now  call  s  and  t  for 
"singlet"  and  "triplet,"  respectively,  rather  than  X  and  u)  conse¬ 
quently,  there  may  be  a  vibrational  overlap  between  adjacent  vibra¬ 
tional  levels  in  N20  which  correspond  asymptotically  to  adjacent 
vibrational  levels  of  N2» 

Using  first-order  perturbation  theory  for  the  vibrational  wave 
function  0^  (g;  R)  of  Eq.  (7), 


\i<£SE>  =  9xi(£)  + 


k  i  i 


exk<£> 


(16) 


13 


Where 


su<£> 


.Urn 

R-*» 


eu<£i  £> 


is  tiie  wave  function  for  vibrational  quantum  number  of  Nji  evi¬ 
dently  there  is  no  dependence  of  X(=  s,t)  in  the  asymptotic  limit. 

v*  .  v( singlet)  Qr  y( triplet) 

( s  or  t ) 

represents  thfi  collisional  perturbation  of  O'  7  on  N2,  as  a 
function  of  j)  and  R. 

Inserting  the  perturbation  ansatz  of  Eq.  (16)  in  the  definition, 
3q.  (12b),  of  qC , 

J^st;  ij(R*j  ~  f  Qsite*  &  etj^J  «*-/«*  etj^d£  + 

I 

<£*£  +  5  ^4^  f  esi(£)8tt(£W£-  (»> 

* 

If  the  collisional  perturbation  V  is  assumed  to  be  linear  in  the 
vibrational  coordinate  p,  the  relations 


Vs  *  ag(p  -  Pe)Us(R) 


V*  *  as(p  -  pe)  Ufc(R) 


(18a) 

(18b) 


follow  for  V  ,  and  to  lowest  order 

1/2 


q  atj  ij 


,00 


0  St,  ii<R> 

1/2 

(i*j) 


*  j  9i9i 


d£=  1 


at  Ut(R)  .  «s  US(R) 

— ErrTj — 


(I|p|i) 


(19a) 


(19b) 
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If  it  is  further  assumed  that  0^(£)  can  be  described  by  harmonic 
oscillator  wave  functions  so  that  only  adjacent  levels  couple, 
(i|p|i-l)  =  /ih/2mu  ,  where  ^  =  reduced  mass  of  N2  and  uj  =  vibra¬ 
tional  (angular)  frequency,  so  that 


U5(R) 


>/i/2gftu)3 


f  20) 


The  harmonic  oscillator  matrix  element  can  be  easily  evaluated 
for  Nj  to  give 

Vi/2phiu3  =  /I/(2.05  x  10“ 3)  cm/erg,  (21) 


while  the  triplet  parameters  can  be  estimated  from  the  results  of 
O-N  scattering  (Ref.  24)  to  give  at  =  2.65A-1  and  for  an  0-N  distance 
of  1.8&  (See  Fig.  1) ,  ll*"  ss  1.5  ev  so  that 

at;Ut(R  1.8  A)  =6.4  x  10“4  erg/cm  .  (22) 

If  it  is  reasonable  to  assume  that  the  singlet  interaction  is  some 
fraction  of  the  triplet,  say  K,  then 

r  p  I1/2 

hstsl,  0<R)  =<1-K)(0.31)  (23) 

For  the  numerical  calculations  that  follow,  K  is  varied  from  0.8  to 
0.2.  As  a  reference  for  discussion,  K  =  0.6  is  taken  as  being  repre¬ 
sentative,  til  us 

1/2 

«  0.124  (24) 

Clearly,  improvements  in  this  number  may  be  made,  but  Eq.  (24)  is 
expected  to  be  roughly  representative  for  the  qualitative  treatment 
which  follows.  The  results  for  other  values  of  K  are  given  in 
Appendix  A. 


c 


st;  1,  0 


(R) 


9> 
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IV 


UNIMOLECULAR  DECOMPOSITION  OF  NjO 


The  nonadiabatic  thermal  decomposition  of  NjO  (*E)  has  been  the 
subject  of  many  investigations  (Refs.  25,  26,  and  27).  Two  experi¬ 
mental  facts  that  emerge  from  these  investigations,  in  the  high- 
pressure  linit,  relate  to  this  curve-crossing  model.  The  first  fact 
is  the  activation  energy  associated  with  the  decomposition  reaction 
of  Eq.  (3), 


N20(xV)  -•  N^X1!*)  ♦  0(3P)  (3) 

thlch  is  found  to  be  60  kcal/mole  (~  2.6  ev).  The  second  is  the 
fact  that  the  pre- exponential  factor  in  the  high-pressure  rate  co¬ 
efficient  i3  from  one  to  two  orders  of  magnitude  smaller  than  the 
expected  gas  kinetic  value.  Both  of  these  results  can  be  qualitatively 
explained  through  the  curve-crossing  model  outlined  in  the  previous 
sections.  The  measured  activation  energy  gives  insight  into  the 
position  of  the  singlet  and  triplet  crossing  point— representing  the 
most  probable  exit  channel  for  the  decomposition.  It  is  not  entirely 
clear  whether  the  decomposition  reaction  occurs  primarily  through  a 
stretching  vibration,  giving  insight  into  the  collinear  configuration 
case,  or  whether  the  decomposition  involves  a  combination  of  stretch¬ 
ing  and  bending  modes.  The  latter  case  might  be  more  important  be¬ 
cause  of  its  larger  statistical  weight;  thus,  the  measured  activation 
energy  presumably  provides  some  angle-averaged  minimum  singlet- 
triplet  crossing  energy. 

The  low  pre-exponential  factors  can  be  explained  by  the  spin 
change  required  between  products  and  reactants,  i.e.,  a  change  from 
the  singlet  reactant  configuration  to  the  triplet  product  configura¬ 
tion  via  a  direct  curve  crossing  (Ref.  26).  On  the  assunption  that 

16 


this  crossing  can  be  characterized  by  the  Landau-Zener  formula,  the 
spin-orbit  interaction  necessary  to  explain  the  low  pre-exponential 
factor  is  of  the  order  of  200  cm”^.  This  value  is  roughly  consistent 
with  the  value  necessary  to  correlate  the  measured  quenching  rate 
coefficient  data  on  OC^D)  by  N0  and  also  to  explain  the  efficiency  of 
0(  P)  in  relaxing  vibrational  energy  in  Ng*  These  subjects  are  the 
content  of  tdie  next  two  sections. 
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V.  QUENCHING  OF  O^)  BY  N2 

Having  established  that  a  direct  curve  crossing  exists  between 
the  singlet  reactant  curves  and  the  triplet  product  curves  as  shown 
in  Figs.  1  and  2,  an  estimate  is  now  developed  of  the  magnitude  of 
the  quenching  cross  section  and  the  likelihood  of  vibrational  energy 
being  channeled  into  the  product  N2<  To  accomplish  this  estimate, 
the  potential  curves  as  shown  in  Figs.  1  and  2  are  expanded  to  include 
the  vibrational  structure  of  both  reactant  and  product  N2  molecules. 
These  molecular  curves  are  shown  in  Figs.  4  and  5  for  the  collinear 
and  120°  configuration  cases. 

In  order  to  perform  the  quenching  calculation,  it  is  assumed,  as 
in  earlier  calculations  (Ref.  1), that  the  curve- crossing  probability 
at  each  intersection  point  is  characterized  by  the  Landau-Zener  equa¬ 
tion,  i.e.,  Eq.  (3).  The  effective  force  terms  in  the  denominator  of 
Tf- see  Eq.  (10)— have  been  estimated  directly  from  the  slopes  at  the 
curve-crossing  intersection  points  as  shown  in  Figs.  4  and  5.  The 
values  used  are  given  in  Table  1. 

The  collision-modified  vibrational  overlap  at  each  crossing  point, 
as  outlined  in  Section  HIC,  has  been  assumed  for  calculational  pur¬ 
poses  to  be[qC]r=  0.124  which  results  when  the  singlet  interaction, 

a  Us,  is  assumed  to  be  60X  of  the  triplet  interaction,  as  shown  in 
s 

Eq.  (24).  The  results  for  other  assumed  values  are  given  in  Appendix  A. 

The  electronic  part  of  the  interaction  matrix  element  at  each 
crossing  point  (the  spin-orbit  interaction)  has  been  treated  as  a 
parameter  of  the  order  of  magnitude  given  for  the  isolated  0(3P)  and 
0(-4)),  as  discussed  in  Section  HIB— see  Eq.  (13b) — and  also  by  the 
results  of  the  unimolecular  decomposition  measurements. 
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FIGURE  5.  Schematic  Potential  Energy  Curves  for  the  N^-O  System,  for  Vibrational- 

A  Q 

Electronic  Coupling:  Bent  Configuration,  Corresponding  toas  N-N-0  =  120 
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TABLE  1.  EFFECTIVE  FORCE  TERMS  USED  IN  TOE 
QUENCHING  OF  0(1D)  BY  Ng 


Triplet 

Vibrational 

Collinear 

120°  Case 

Level 

Case,  ev/A 

ev/A 

6 

0.7 

0.65 

5 

1.29 

1.03 

4 

1.8 

1.26 

3 

2.0 

1.5 

2 

2.6 

1.54 

1 

2.92 

1.58 

0 

4.28 

1.5 

Following  the  definition  of  the  curve- crossing  probability  at 
each  intersection  point  as  outlined  above,  a  classical  diffusion 
calculation  has  been  performed  as  a  function  of  relative  kinetic 
energy.  The  details  of  this  diffusion  calculation  have  been  pre¬ 
sented  elsewhere  (Refs.  1  and  28). 


Hie  results  for  the  partial  quenching  cross  sections  for  the 
collinear  and  the  120°  configuration  cases  are  shown  in  Tables  2 
a  id  3,  respectively.*  For  these  calculations  the  effective  force 
t  ?ms  listed  in  Table  1  have  been  used.  As  pointed  out  in  Section 
IV,  the  singlet-triplet  crossing  point  was  established  from  the  60 
kcal/mole  activation  energy  determined  in  unimolecular  decomposition 


"ah  of  these  cross  sections  refer  to  the  collision  partners  N2(X^£g) 

and  OC3!))  entering  along  the  lowest  possible  potential  energy  curve, 

i.e.,  the  1£+  curve  of  Fig.  1  or  the  '  curve  of  Fig.  2.  Thus,  the 
total  cross  sections  must  be  multiplied  by  the  statistical  weight 
factor  g  *  0.2 — see  discussion  below  Eq.  (25). 


«*• 
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TABLE  2.  PARTIAL  QUENCHING  CROSS  SECTIONS  FOR 
OC1!))  +  N2  REACTION  (IN  10"16  <2M2)* 

COLLI NEAR  CONFIGURATION 


** 


Level  \vel 

100  cm”1 

300 

500 

6 

0.0 

0.0 

0.0 

5 

0.0 

0.0 

0.0 

4 

0.0 

0.0 

0.0 

3 

0.0 

0.0 

0.0 

2 

0.0 

0.0 

0.0 

1 

0.0 

0.0 

0.12 

0 

0.41 

3.33 

7.57 

Total  Quenching 

Cross  Section 

0.41 

3.33 

7.69 

ft 

Franck  Condon  Factor  =  0.015 
Relative  kinetic  energy  =  0. 

(Corresponds  to  K  = 
1  ev. 

0.6  in  Eq 

Assumed  spin- orbit  interaction. 

TABLE  3.  PARTIAL  QUENCHING  CROSS  SECTIONS  FOR 
0(XD)  +  N2  REACTION  (IN  10"16CM2)* 

BENT  (120°)  CONFIGURATION 

Vibrational  \  v  ** 

Level  \vel 

100  cm”1 

300 

500 

6 

0.0 

0.0 

0.0 

5 

0.0 

0.0 

0.0 

4 

0.0 

0.0 

0.0 

3 

0.0 

0.0 

0.0 

2 

0.0 

0.0 

0.02 

1 

0.0 

0.03 

0.66 

0 

1.54 

9.42 

12.20 

Total  Quenching 

Cross  Section 

1.54 

9.45 

12.88 

w 

Franck  Condon  Factor  =  0.015 

(Corresponds  to  K  = 

0.6  in  Eq. 

Relative  kinetic  energy  *  0.1  ev. 
Assumed  spin-orbit  interaction. 
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measurements.  The  curves  shown  in  Figs.  4  and  5  are  based  on  this  val¬ 
ue,  i.e.,  the  energy  separation  from  the  minimum  well  depth  of  the  N20 
molecule  to  the  crossing  of  the  v  =  0  singlet  and  triplet  curves,  and 
the  effective  force  terms  shown  in  Table  1  have  been  estimated  from 
these  figures.  It  is  unlikely,  however,  that  the  unimolecular  decom¬ 
position  products  are  vibrationally  unexcited  because  the  transition 
probably  passes  through  a  vibrationally  excited  intermediate  state. 
Therefore  the  measured  activation  energy  probably  represents  the 
energy  separation  between  the  minimum  NjO  well  depth  and  some  excited 
vibrational  level  of  the  triplet  curve.  “Die  net  effect  of  this  would 
be  to  reduce  the  effective  force  terms  at  each  crossing  point  from 
those  given  in  Table  1.  In  order  to  test  the  effect  of  this  on  the 
partial  quenching  cross  sections,  the  effective  force  terms  shown  in 
Table  1  have  been  arbitrarily  reduced  by  30%  and  the  <  ■  aching  cal¬ 
culation  repeated  for  both  the  collinear  and  120°  configuration  cases. 
!Hie  results  for  these  calculations  are  given  in  Tables  4  and  5.* 

Several  features  of  these  results  are  of  importance: 

a.  The  total  quenching  cross  section  increases  with  the 
magnitude  of  the  assumed  spin- orbit  interaction. 

b.  The  total  quenching  cross  section  is  larger  for  the 
120°  case  than  for  the  collinear,  as  expected  from 
the  reduced  effective  force  terms  for  the  bent  con¬ 
figuration  geometry. 

c.  The  quenching  cross  sections  are  only  weakly  depend¬ 
ent  on  collision  energy  due  to  the  attractive  nature 
of  the  interaction  potential  for  the  reactants. 

d.  The  interaction  matrix  element  necessary  to  give  rough 
agreement  with  experimental  quenching  measurements  is 

^  | 

All  of  these  cross  sections  refer  to  the  collision  partners  N2(X  Eg) 

and  0(h))  entering  along  the  lowest  possible  potential  energy  curve, 

i.e.,  the  1E+  curve  of  Fig.  1  or  the  1A '  curve  of  Fig.  2.  Thus,  the 
total  cross  sections  must  be  multiplied  by  the  statistical  weight 
factor  g  s  0.2 — see  discussion  below  Eq.  (25). 
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TABLE  4.  PARTIAL  QUENCHING  CROSS  SECTIONS  FOR 


0(1D)  +  N 

REACTION  (IN  10' 

16  CM2)* 

COLLI NEAR  CONFIGURATIONS 

Vibrational  \  *** 

Level  \vel 

100  cm-1 

300 

500 

6 

0.0 

0.0 

0.0 

5 

0.0 

0.0 

0.1 

4 

0.0 

0.0 

0.0 

3 

0.0 

0.0 

0.0 

2 

0.0 

0.0 

0.0 

1 

0.0 

0.0 

0.21 

0 

0.58 

4.54 

9.44 

Total  Quenching 

Cross  Section 

0.58 

4.54 

9.65 

Franck  Condon  Factor  =  0.015  (Corresponds  to  K  =  0.6  in  Eq.  23) 
^Relative  kinetic  energy  *  0.1  ev. 

0.7  times  effective  force  terms  shown  in  Table  1. 

Assumed  spin-orbit  interaction. 


TABLE  5.  PARTIAL  QUENCHING  CROSS  SECTIONS  FOR 
OC^)  +  N.  REACTION  (IN  10“16  cm2)* 


BENT* 

£ 

(120°)  CONFIGURATION 

Vibrational 

Level 

\  *** 

\Vel  100  cm-1 

300 

500 

6 

0.0 

0.0 

0.0 

5 

0.0 

0.0 

0.0 

4 

0.0 

0.0 

0.0 

3 

c.c 

0.0 

0.0 

2 

0.0 

0.01 

0.03 

1 

0.0 

0.23 

1.22 

0 

2.16 

11.14 

10.64 

Total  Quenching 
Cross  Section 

2.16 

11.38 

11.89 

* 

Franck  Condon  Factor  = 

0.015  (Corresponds  to  K  = 

0.6  in  Eq 

** 


Relative  kinetic  energy  =  0.1  ev. 

0.7  times  the  effective  force  listed  in  Table  1. 
Assumed  spin-orbit  interaction. 
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consistent  with  the  unimolecular  decomposition  calculation 
presented  in  the  previous  section,  but  this  value  is  sig¬ 
nificantly  larger  than  the  spin-orbit  interaction  energy 
for  an  isolated  oxygen  atom  (cf.  Ref.  21). 
e.  Due  to  the  small  vibrational  overlap  between  adjacent  vi¬ 
brational  levels,  the  partial  cross  sections  for  exciting 
N2  vibrational  levels,  v  >0,  is  small;  and  the  fraction 
of  collisions  leading  to  the  v  >  1  exit  channel  is  of  the 
order  of  0.10  or  smaller  (note  that  the  fraction  of  energy 
which  channels  into  excited  vibrational  levels  of  N2  is 
roughly  (0.1) (1/7)  <  2%). 


Here  the  effective  de-excitation  cross  section,  ad>  has  been  es 
timated  and  found  to  be  only  weakly  dependent  on  velocity  u  (or  tem¬ 
perature  T).  The  expression  for  the  de-excitation  rate  constant 
kd(T)  is 


kd(T)  =  g 


J ad(u)  f(u,T)  du, 


(25a) 


where  f(u,  T)  is  the  (Maxwell- Boltzmann)  velocity  distribution  func¬ 
tion  and  g  is  a  statistical  weight  factor;  because  ad  is  only  weakly 
dependent  on  u,  kd(T)  may  be  reduced  to  the  following: 

kd(T)  =  ST  (25b) 

vdiere  u(T)  is  the  thermal  mean  speed  8  x  104  cm/sec  at  T  =  300°K. 

The  statistical  weight  factor,  g,  is  the  probability  that  the 
0(b)  atom  approaches  the  N2  molecule  in  the  configuration  corresponding 
to  a  crossing.  For  the  collinear  case  (cf .  Fig.  1),  this  corresponds 
to  the  state  rather  than  the  or  ^n»  i.e..  g  =  1/(1  42+2)=  1/5. 
For  noncollinear  configurations,  the  statistical  factor  is  uncertain, 
since  additional  curve  crossings  may  take  place,  involving  the  two 
repulsive  states.  Thus,  using  g  =  0.2  (which  may  be  an  underestimate) 
and  u(T)  =  8  x  10  4  cm/sec,  the  range  of  experimental  values  of  kd(T), 
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which  is  (4-9)  x  10-11  cm3/sec  (Ref.  13),  corresponds  to  tj.  =  (2.5  - 

-15  2  a 

5.6)  x  10  cm  ,  or  two  to  six  times  the  calculated  value  of  the 

cross  section.  In  fact,  the  lower  experimental  value  is  more  plausi¬ 
ble,  and  agreement  with  experiment  within  a  factor  2-3  is  certainly 
all  that  can  be  expected  by  a  calculation  as  primitive  as  the  present 
one. 
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VI.  N2  VIBRATIONAL  RELAXATION  BY  0(3P)  atoms 

The  results  of  the  previous  two  sections  can  be  used  to  calculate 
the  vibrational  relaxation  cross  sections  for  N2  by  0  atoms  by  analogy 
to  the  curve- crossing  model  for  the  relaxation  of  N2  by  alkali  atoms 
(Ref.  3).  For  an  estimate  of  this  process,  refer  to  Fig.  4,  the 
collinear  interaction  of  0(3P)  atoms  with  N^1^*),  i.e.,  reaction  (5). 

As  can  be  seen  from  this  figure,  this  mechanism  for  atom  relaxa¬ 
tion  is  accompanied  by  an  activation  energy  of  about  0.8  ev  from  the 
reactants  to  the  first  crossing  point  of  the  v  =  0  singlet  surface. 

The  activation  energy  associated  with  this  collinear  crossing  con¬ 
figuration  will  be  considerably  lower  than  for  noncollinear  config¬ 
urations  due  to  the  decrease  in  the  well  depth  of  the  NgO  singlet 
surface  for  angles  less  than  180°.  It  should  be  remarked,  however, 
that  since  the  unimolecular  decomposition  of  NjO  may  preferentially 
give  vibrationally  excited  N2,  an  activation  energy  of  0.8  ev  is 
likely  to  be  a  maximum  value. 

The  probability  of  transition  from  the  v  =  1  reactant  surface  to 
the  v  =  0  product  surface,  P,  is  roughly  given  by  2P  (1  -  P.  )(1  -  P  ), 
vhere  Pfl  and  P^  are  the  diabatic  (i.e.,  curve-crossing)  probabilities 
at  the  v  =  0,  v  =  0  and  v  =  1,  v  =0,  intersections,  respectively. 
The  probabilities  will  be  considerably  larger  than  in  the  corre¬ 
sponding  quenching  case  since  the  velocity  at  the  crossing  points 
will  be  small  by  comparison.  Very  crudely,  the  vibrational  deactiva¬ 
tion  probability  may  be  taken  as  the  ratio  of  the  v  =  1  to  v  =  0 
partial  quenching  cross  sections  since  the  increased  probabilities 
due  to  a  considerably  lower  velocity  at  the  crossing  point  will  be 
counteracted  by  the  additional  adiabatic  curve  crossing  necessary  to 
transfer  from  the  v  =  1  reactant  surface  to  the  v  =  0  singlet  surface. 
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For  the  assumed  range  of  spin-orbit  interaction  values,  the  ratio  of 
the  v  =  1  to  v  =  0  cross  sections  is  roughly  0.1.  The  relaxation 
rate  coefficient  is  thus 

“Eact^T 

p10  =  spe  »  (26) 


vhere 

steric  factor  to  account  for  noncollinear  configurations 
1/6 

collisional  deactivation  probability/collision 

?.(v-= n..  -  o.i 

o(v  =  0) 


activation  energy  to  reach  first  curve  crossing 
0.8  ev. 

This  numerical  value  for  the  collisional  deactivation  probability 
per  collision  is  confirmed  by  more  detailed  calculations  but  is 
clearly  sensitive  to  the  vibrational  overlap  and  the  assumed  spin- 
orbit  interaction  at  the  crossing  points,  as  in  the  quenching  case. 

The  important  point  is  that  the  temperature  dependence  is  weak  in 
the  pre- exponential  factor  and  should  mainly  come  from  the  activation 
energy  term.  Inserting  the  above  numbers  into  Eq.  (26),  we  get 

P1Q  =  1.7  x  10“ 2  exp  £-9.5  x  103/T(°K))  .  (27) 

In  order  to  compare  this  value  with  the  available  experimental 
data,  the  probability  can  be  converted  to  a  relaxation  time,  t,  by 


-hcu)/kT ' 

>  > 


(28) 
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vfaere  Z  —  collision  frequency 


uu  =  vibrational  spacing 
Nfi  =  number  density 

and  t  =  relaxation  time 

Using  a  =  3.02  A,  we  get 

^t4000°K  =  ^  sec  atm 

and 

^Tjqqqoj^  1*1  |jl  sec  atm 

where  P  is. the  pressure  (in  atmospheres). 

Figure  6  shows  the  calculated  relaxation  time  for  the  relaxation 
of  N2  vibrational  energy  in  collisions  with  0  atoms  as  compared  with 
available  data.  For  extrapolation  to  other  temperatures,  the  equation 
of  the  calculated  relaxation  time  is  (T  in  °K,  P  in  atm,  T  in  sec) 


(1/Pt)  =  1.64  x  109  T-1/2  e-950°/T  [1  -  e-3250/1] 


(29) 


It  must  be  stressed  that  the  "non-adiabatic"  or  ’’chemical"  mecha¬ 
nism  for  the  vibrational  relaxation  of  N2  by  0  that  has  been  discussed 
here  is  not  the  only  one.  There  is  also  the  ’’adiabatic"  or  "Landau- 
Teller"  mechanism — cf .  Ref.  29.  Using  the  semi- empirical  model  of 
Millikan  and  White  (kef.  30)  for  the  Landau- Teller  mechanism,  we  find 
that 


(Pt; 


chem 
(see  Eq. 


29) 


<PT>Landau-Teller  for  T  2  600  K 


9- 


29 


(30) 


4500  4000  3500  3000  2700 
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FIGURE  6.  Vibrational  Relaxation  of  Nj  by  Collision  with  O.  The  Theoretical 
Curve  of  Eq.(29)  is  Superposed  on  the  Experimental  Data  of 
Breshears  and  Bird  (Reft  14). 
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In  other  words,  for  T  >  600°K  the  chemical  mechanism  dominates,  while 
for  T  <  600°K  the  0.8  ev  activation  energy  for  the  chemical  mechanism 
presents  such  a  big  barrier  that  the  nonchemical  Landau-Teller  mechan¬ 
ism  dominates.  In  fact,  a  comparable  situation  is  known  to  occur  for 
NO-NO  collisions— see  Refs.  22  and  31. 
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VII.  DISCUSSION 


From  the  results  of  the  previous  sections,  several  features  of 
the  various  processes--in  particular  of  the  quenching  reaction  (4)— 
can  be  indicated. 

1.  It  may  seem  surprising  that  the  general  order  of  magnitude 
of  these  rather  large  cross  sections  [*.  (1/10)  gas  kinetic]  can  be 

pi  _■] 

explained  by  the  small  interaction  V  ~  100-500  cm  (0.01  -  0.04  ev) 
associated  with  the  spin-change  transition.  This  follows  immediately 
from  the  assumption  of  a  Landau-Zener  transition,  i.e.,  on  interactions 
occurring  over  relatively  small  distances  and  calculable  by  perturba¬ 
tion  theory.  There  is  the  possibility  that  these  assumptions  are  in¬ 
adequate,  so  that  "close-coupling"  effects  may  be  important.  It 
should  be  noted  that  close-coupling  calculations  for  a  problem  as  com¬ 
plex  as  the  present  one  are  hardly  feasible  at  present;  nor  are  they 
likely  to  become  feasible  in  less  than  three  to  five  years. 

2.  The  vibration-electronic  coupling  in  this  reaction  appears 
to  be  very  small.  This  result  is  a  direct  consequence  of  the  small 
overlap  between  adjacent  vibrational  levels  even  in  the  collisionally 
perturbed  case.  Although  the  estimate  of  the  modified  Franck-Condon 
factor  is  crude,  it  is  rather  unlikely  that,  even  with  a  more  adequate 
estimate  of  the  vibrational  overlap,  a  significant  overlap  between 
vibrational  levels  separated  by  more  than  one  quantum  will  be  pre¬ 
dicted.  This  point  is  made  even  stronger  by  the  relatively  large 
cross  sections  predicted  for  the  v  =  0  channel,  i.e.,  no  vibrational 
excitation.  It  is  anticipated  that  other  quenching  reactions  which 
proceed  entirely  via  the  crossing  of  covalent  potential  energy  curves 
will  show  similar  behavior.  This  conclusion  is  in  marked  contrast  to 
the  results  of  quenching  via  ionic  intermediate  complexes  where  the 
vibration-electronic  coupling  is  large. 
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3.  In  the  calculations  presented  here  only  one  of  the  triplet 
curves  (see  Fig.  1)  was  included  explicitly.  Inclusion  of  the  second 
curve  could  increase  the  interaction  matrix  element  for  crossing 
points  greater  than  the  v  =  0  level;  however,  there  appears  to  be  no 
adequate  way  of  accounting  for  the  interaction  at  the  intersection  of 
three  potential  energy  curves . 

In  addition,  the  position  of  the  crossing  point  of  the  v  =  0 
singlet  and  triplet  curves  was  established  by  the  activation  energy 
of  the  high-pressure  unimolecular  decomposition  of  N20,  as  discussed 
in  Section  IV.  It  is  likely  that  this  activation  energy  would  be 
more  representative  of  a  higher  lying  vibrational  level,  e.g.,  the 
intersection  of  the  v  =  3  singlet  and  triplet  curves,  and  thus  the 
actual  effective  force  terms  given  in  Table  1  would  be  reduced  with 
a  corresponding  increase  in  the  adiabatic  curve-crossing  character 
and  hence,  total  quenching  cross  sections. 

This  conclusion  is  verified  by  the  calculations  based  upon  0.7 
times  the  effective  force  terms  as  given  in  Table  1. 

4.  The  maximum  contribution  to  the  total  quenching  cross  section 
appears  to  come  from  noncollinear  configuration  trajectories.  This 
conclusion  is  a  result  of  the  reduction  in  the  effective  force  evident 
at  the  curve-crossing  points  due  to  the  decrease  in  the  singlet  well 
depth  with  decreasing  angle  a  =  N  -  N  -  0  from  180°.  As  a  consequence 
of  this  reduction,  the  N2  quenching  product  is  likely  to  be  rotationally 
excited  following  a  successful  quenching  reaction. 

The  degree  of  rotational  excitation  may  be  estimated  as  follows. 

For  an  N-N  internuclear  distance  rg,  the  effective  impact  parameter 
b  of  the  collision  is 


b  =  %  rg  sin  (rr  -  a) 


(30a) 


For  incident  velocity  u,  the  initial  angular  momentum  available  for 

transfer  to  rotation  is  Mn>w  ub  corresponding  to  a  maximum  rotational 

u  ?  W2 

quantum  number 
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(30b) 


max 


=  M 


0;N 


ub/h 


where 


M0;N,,  M0  MN2^M0  +  MN2^  ~  1,69  X  10  23  gm* 

Using  u  =  4.5  x  10  9  cm/sec  (corresponding  to  a  kinetic  energy  of  1  ev, 
which  is  appropriate  at  the  crossing  point— see  Figs.  1  and  2)  and 
b  =  0.5  x  10'8  cm  gives  Jmax  =  35,  so  that  max  =  Be  J(J+1)  0.31  ev; 

thus,  at  most  15  percent  of  the  excitation  energy  (=  2  ev)  can  go  into 
rotation. 

In  s’lmmary,  on  the  basis  of  qualitative  features  of  the  potential 
energy  surfaces  for  the  quenching  of  0(^D)  by  N2,  less  than  5  percent 
of  the  electronic  quantum  of  O('h))  is  expected  to  channel  into  the 
vibrational  levels  of  N2.  Thus,  reaction  (2)  will  not  be  a  significant 
source  of  vibrational  energy  in  the  upper  atmosphere.  In  addition,  the 
relaxation  of  N2  vibrational  energy  by  0(3P)  atoms  at  high  temperatures 
is  consistent  with  this  curve-crossing  model  and  suggests  a  strong  tem¬ 
perature  dependence  for  extrapolation  to  lower  temperatures  appropriate 
to  upper  atmospheric  conditions.  Note,  however,  that  below  about  600°K 
the  effect  of  "adiabatic"  or  "Landau- Teller"  vibrational  relaxation  of 
N2  by  0  is  likely  to  predominate  over  the  present  "chemical"  mechanism. 
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APPENDIX  A 


PARTIAL  QUENCHING  CROSS  SECTIONS  FOR 

1  o(v)  7 

0(1))  +  N2  — ►  0(  P)  +  N2(v) 

AS  A  FUNCTION  OF  THE  COLLISION-MODIFIED  FRANCK "CONDON  FACTOR* 


All  of  the  cross  sections  in  this  Appendix  refer  to  thr  collision 
partners  N2(X12^)  and  0(*D)  entering  along  the  lowest  possible  poten¬ 
tial  energy  curve,  i.e.,  the  ^E+  curve  of  Fig.  1  or  the  ^A '  curve  of 
Fig.  2.  Thus,  the  total  cross  sections  must  be  multiplied  by  the 
statistical  weight  factor  g  >  0.2 — see  discussion  below  Eg.  (25). 
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TABLE  A-l.  PARTIAL  QUENCHING  CROSS  SECTIONS  FOR 
O^D)  +  N2  REACTION  (IN  10“ 16  CM2)* 

BENT  (120°)  CONFIGURATION 


Vibrational  \  ..  ** 
Level  \vel ’ 

100  cm-1 

300 

500 

6 

0.0 

0.0 

0.0 

5 

0.0 

0.0 

0.0 

4 

0.0 

0.0 

0.0 

3 

0.0 

0.0 

0.0 

2 

0.0 

0.0 

0.0 

1 

0.0 

0.0 

0.07 

0 

1.54 

9.50 

12.56 

Total  Quenching 

Cross  Section 

1.54 

9.50 

12.63 

Franck  Condon  Factor  =  0.004  (Corresponds  to  K  =  0.8  in  Eq.  23). 
Relative  kinetic  energy  =0.1  ev. 

** 

Assumed  spin-orbit  interaction. 


TABLE  A- 2.  PARTIAL  QUENCHING  CROSS  SECTIONS  FOR 

.1  .  _  .  R  9  * 


0(1D)  +  N_  REACTION  (IN  10  16  CM2) 

2  n  ** 

BENT  (120°)  CONFIGURATION 

Vibrational  \  *** 

Level  \  vel 

100  cm-1 

300 

500 

6 

0.0 

0.0 

0.0 

5 

0.0 

0.0 

0.0 

4 

0.0 

0.0 

0.0 

3 

0.0 

0.0 

0.0 

2 

0.0 

0.0 

0.0 

1 

0.0 

0.0 

0 .20 

0 

2.16 

11.33 

11.05 

Total  Quenching 

Cross  Section 

2.16 

11.33 

11.25 

* 

Franck  Condon  Factor  =  0.004  (Corresponds 
Relative  kinetic  energy  =  0.1  ev. 

** 

0.7  times  effective  force  terms  shown  in  ' 

to  K  =  0 

Table  1. 

.8  in  Eq.  23) 

For  120° 

configuration  case. 

*** 

Assumed  spin- orbit  interaction. 
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TABLE  A- 3.  PARTIAL  QUENCHING  CROSS  SECTIONS  FOR 
Oib)  +  N,  REACTION  (IN  10-16  CM2)* 

COLLINEAR 

CONFIGURATION 

VibrationalX  ** 

Level  \vel 

100  cm”1 

300 

500 

6 

0.0 

0.0 

0.0 

5 

0.0 

0.0 

0.0 

4 

0.0 

0.0 

0.0 

3 

0.0 

0.0 

0.0 

2 

0.0 

0.0 

0.0 

1 

0.0 

0.0 

0.0 

0 

0.41 

3.33 

7.60 

Total  Quenching 

Cross  Section 

0.41 

3.33 

7.60 

* 

Franck  Condon  Factor  =  0.004 

(Corresponds 

to  K  =  0.8 

in  Eq. 

Relative  kinetic  energy  =  0.1  ev. 
Assumed  spin-orbit  interaction. 


TAB.E  A-4.  PARTIAL  CLENCHING  CROSS  SECTIONS  FOR 
O^D)  +  N„  REACTION  (10"16  CM2)* 


\ 


el 


Vibrational 
Level 

6 

5 

4 

3 

2 

1 

0 

Total  Quenching 
Cross  Section 


COLLI NEAR  CONFIGURATION 

*** 


)0  cm”1 

300 

500 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.58 

4.54 

9.52 

0.58 

4.54 

9.52 

Franck  Condon  Factor  =  0.004  (Corresponds  to  X  =  0.8  in  Eq.  23) 
Relative  kinetic  energy  -  0.1  ev. 

0."  Times  effective  force  terms  shown  in  Table  1  for  collinear 
configuration  case. 

Assumed  spin-orbit  interaction. 
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TABLE  A- 5.  PARTIAL  QUENCHING  CROSS  SECTIONS  FOR 
O^D)  +  n2  REACTION  (10“16  CM2)* 

BENT  (120°)  CONFIGURATION 


VibrationaA 


Level 

\  el 

100  cm”1 

300 

500 

6 

0.0 

0.0 

0.0 

5 

0.0 

0.0 

0.0 

4 

0.0 

0.0 

0.0 

3 

0.0 

0.0 

0.0 

2 

0.0 

0.01 

0.08 

1 

0.0 

0.33 

1.37 

0 

1.54 

9.25 

11.76 

Total  Quenching 

Cre  ss  Section 

1.54 

9.59 

13.21 

Franck  Condon  Factor  =  0.035  (Corresponds  to  K  =  0.4  in  Eq.  23). 
Relative  kinetic  energy  -  0.1  ev. 

Assumed  spin-orbit  interaction. 


TABLE  A-6 .  PARTIAL  QUENCHING  CROSS  SECTIONS  FOR 
CO(1D )  +  N„  REACTION  (10'16  CM2)* 

BENT  (120;  CONFIGURATION 


Vibrational\  *** 
Level  \vel 

6 

5 

4 

3 

2 

1 

0 

Total  Quenching 
Cross  Section 


10  cm”1 

300 

500 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.03 

0.12 

0.0 

0.50 

2.46 

2.16 

10.91 

10.11 

2.16 

11.44 

12.69 

Franck  Condon  Factor  =  0.035  (Corresponds  to  K  =  0.4  in  Eq.  23). 
Relative  kinetic  energy  =  0.1  ev. 

0.7  times  effective  force  terms  shown  in  Table  1  for  120c 
configuration  case. 


Assumed  spin-orbit  interaction. 
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TABLE  A-7.  PARTIAL  QUENCHING  CROSS  SECTIONS  FOR 

_  ,1_  .  „  _ _ -IK  0  * 


0(1D)  +  N 

REACTION  (10' 

16cmV 

COLLINEAR  CONFIGURATION 

Vibrational\.;  ** 
Level  \  el 

100  cm-1 

300 

500 

6 

0.0 

0.0 

0.0 

5 

0.0 

0.0 

0.0 

4 

0.0 

0.0 

0.0 

3 

0.0 

0.0 

0.0 

2 

0.0 

0.0 

0.01 

1 

0.0 

0.0 

0.38 

0 

0.41 

3.33 

7.46 

Total  Quenching 

0.41 

3.33 

7.85 

Cross  Section 

Franck  Condon  Factor  =  0.035  (Corresponds  to  K  =  0.4  in  Eq.  23). 
Relative  kinetic  energy  =0.1  ev. 

#* 

Assumed  spin-orbit  interaction. 


TABLE  A- 8. 


PARTIAL  QUENCHING  CROSS  SECTIONS  FOR 
O^D)  +  N_  REACTION  (10"16  CM2)* 
COLLINEAR  CONFIGURATION 


Vibrational \  *** 

Level  \vel 

6 

5 

4 

3 

2 

1 

0 

Total  Quenching 
Cross  Section 


10  cm"1 

300 

500 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.02 

0.0 

0.22 

0.50 

0.58 

4.51 

9.26 

0.58 

4.73 

9.78 

* 

Franck  Condon  Factor  =  0.035  (Corresponds  to  K  =  0.4  in  Eq.  23). 
Relative  kinetic  energy  =  0.1  ev. 

0.7  times  effective  force  terms  shown  in  Table  A-l  for  collinear 

configuration  case. 

nUr 

Assumed  spin-orbit  interaction. 
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